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The quintessential property of neuronal systems is their intensive patterns of selective synaptic 
connections. The current work describes a physics-based approach to neuronal shape modeling 
and synthesis and its consideration for the simulation of neuronal development and the formation 
of neuronal communities. Starting from images of real neurons, geometrical measurements are 
obtained and used to construct probabilistic models which can be subsequently sampled in order to 
produce morphologically realistic neuronal cells. Such cells are progressively grown while monitoring 
their connections along time, which are analysed in terms of percolation concepts. However, unlike 
traditional percolation, the critical point is verified along the growth stages, not the density of cells, 
which remains constant throughout the neuronal growth dynamics. It is shown, through simulations, 
that growing beta cells tend to reach percolation sooner than the alpha counterparts with the same 
diameter. Also, the percolation becomes more abrupt for higher densities of cells, being markedly 
sharper for the beta cells. 

PACS numbers: 89.75.Fb, 87.18. Sn, 02.10.Ox, 89. 75. Da, 89. 75. He, 2.50.Ey, 2.50.Ng 



The brain is a world consisting of a number of 
unexplored continents and great stretches 
of unknown territory. (Santiago Ramon- 
y-Cajal) 

Neurons can be understood as cells which, along the 
evolutionary process, have become highly specialized for 
establishing connections between themselves along a wide 
range of spatial scales (ranging from microns to meters). 
In order to minimize metabolism and allow connections 
to selective targets, neurons acquired their intricated, 
ramified shapes. Indeed, instead of implementing ca- 
sual connections with every surrounding cell, a neuron 
links to specific targets which can be nearby in the same 
neuronal region or far away in another cortical area origi- 
nating, in the process, the basic architecture required for 
proper operation of the central nervous system. Inter- 
estingly, the connectivity pattern of a mature neuronal 
system is determined not only by the genome, which is 
unable to code all connections |2J, but predominantly 
by the history of neuronal activity under stimuli presen- 
tation. Neurons are produced at ventricular zones of the 
neuroepithelium, in the form of neuroblasts, which there- 
fore differentiate and migrate to specific target regions 
and start to unfold their dendritic and axonal processes 
0. As such structures develop and extend towards spe- 
cific targets, which occurs under the guidance of trophic 
factors, they synapse and start forming communities (or 
clusters) of connected cells, organized in specific ways so 
as to achieve proper operation. Indeed, the functional 
properties of such structures are to a large extent related 



to the underlying connecting patterns, implying that one 
of the fundamental problems in neuroscience is to under- 
stand how neuronal connections are established during 
the development of the central nervous system pj . 

Since connectivity is the main purpose underlying neu- 
ronal growth and organization, it is interesting to obtain 
suitable mathematical structures and relationships capa- 
ble of representing and modeling the development of neu- 
ronal systems at a high level of morphological realism. 
While graphs/networks, where neurons are assigned to 
nodes and synapses to edges, provide a natural means to 
express the neuronal connections, several concepts from 
statistical mechanics can be used to model and simulate 
the connection dynamics. By providing an interesting in- 
terface between graph theory and statistical mechanics, 
the recent area of complex networks 0- C3- Cfll represents a 
particularly promising perspective to bridge the gap be- 
tween the morphology and dynamics of neuronal systems. 
In particular, the concept of percolation [f| stand out as 
particularly relevant for such investigations. Previous re- 
lated works include the statistical physics investigation of 
scaling properties and the degree of separation in cortical 
networks [g , the small- world characterization of neuronal 
structures grown in vitro the use of critical percola- 
tion point for neuronal shape characterization [8j|, and 
the identification of electrically active clusters in neu- 
ral networks 0. Related works addressing the relation- 
ship between neuronal geometry and function can also be 
found in the literature (e.g. 0, 0, 0,0). While 
such works have considered stactic neuronal shapes, the 




development of a framework to model neuromorphically 
realistic neurons reported in 0] allows investigations of 
the neuronal connectivity during simulated neuronal de- 
velopment, by monitoring the size and other properties 
of the existing clusters in terms of time. Such a perspec- 
tive motivated the extension of the concept of percola- 
tion to consider growing structures where the shapes of 
the objects may vary with time, a possibility proposed 
and investigated possibly for the first time in the present 
work. 

This article starts by describing how the neuronal cells 
are represented and statistically modeled in terms of 
probabilities and follows by presenting the simulation of 
neuronal growth by using the Monte Carlo approach, as 
well as the characterization of the obtained structures in 
terms of the maximum cluster size observed along time. 
Such issues are illustrated with respect to a database of 
2D neuronal cells including cat retina ganglionar cells of 
the types alpha (23 samples) and beta (27 samples), of 
which typical cells are illustrated in Figure 



I. NEURONAL MODELING 

One first key isse in neuromorphic modeling regards 
how to represent mathematically the geometry of neu- 
rons. While the typically observed diversity of shapes 
for the same class of cells immediately implies the use 
of statistics, the choice of the best (in the sense of be- 
ing the most compact) set of morphometric parameters 
capable of representing the neuronal shape without con- 
siderable loss of information remains a challenging issue 
[l3| . The methodology for 2D neuronal representation 
adopted in this work follows the framework reported in 



|14j , involving a probabilistic model considering the num- 
ber of branches, the angles between them, the length of 
the dendritic and axonal segments, branching probability, 
and the length and angle of arcs of each branch. There- 
fore, the first step is to obtain such measurements from 
images of the real neuronal cells to be modeled. Typ- 
ically, the cells are histologically marked and prepared, 
mounted on slides, and the respective images acquired 
through a camera interfaced to a light transmission (or 
fluorescence) microscope. The neurons in such images are 
then identified and isolated (e.g. jl7j|'). producing binary 
representations (i.e. images containing only the neuronal 
cell - marked as one, and the background - marked as 
zero). An alternative way to obtain the binary images 
of the neural cells is through camera-lucida drawings, as 
is the case for the images in Figure ^ Once such binary 
images are obtained, their boundaries are extracted by 
using a customized neural tracer |22j, producing results 
such as those illustrated in Figure which corresponds 
to the boundaries of the cells in Figure ^ 

The dendrites are henceforth understood as trees, so 
that the respective hierarchical level can be precisely de- 
fined while considering the soma as reference. Therefore, 
the dendritic segments directly connected to the soma, as 
well as branching points initiating at such segments, are 
identified as being at hierarchical level 1, and so on. Our 
simulations are restricted to a maximum of 10 hierarchi- 
cal levels, as there are very few branchings occurring at 
higher levels in the real cells. The probability of branch 
points for each considered type of cell are shown in Fig- 
ure l^a) and (b), respectively for alpha and beta cells. 
The probability of the number of dendritic segments di- 
rectly attached to the soma (i.e. hierarchical level 1) is 
also necessary for the statistical model of the growing 
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FIG. 3: Probability function illustrating the branching prob- 
ability of alpha (a) and beta (b) cells. 



cells. Figure 0| show the cumulative densities for the al- 
pha (a) and beta (b) neuronal types. 

In addition to being essential for neuronal shape mod- 
eling, the above branching and initial densities provide 
interesting information by themselves. For instance, it 
is clear from the two densities in Figure [21 that the al- 
pha cells are characterized by higher branching rates at 
the lower hierarchical levels, as is clear from the more 
accentuated decreasal of the respective density along the 
hierarchical levels. Figure 0] presents the probability of 
the number of primary dendritic segments for alpha (a) 
and beta (b) cells. These graphs indicate that the beta 
cells are more likely to branch than alpha cells at hierar- 
chical levels between 8 and 12. 

In addition to the above probabilities, it is also nec- 
essary to obtain probabilistic models of the dendritic 
segment arc-lengths. Although alpha cells are typically 
much larger than beta cells (especially in the periph- 
ery of the retina), we used size- normalized versions of 
the considered neurons in order to have neurons with 
similar sizes. This allows our percolation study to be 
mostly defined by the shape intrincacy of the cells rather 
than their sizes |2a |. Figure [5] shows the cumulative two- 
variated density of such lengths in terms of the hierar- 
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FIG. 4: Probability function of the number of primary den- 
dritic segments to the alpha (a) and beta (b) cells. 



chical level for the alpha (a) and beta (b) types of cells, 
while Figure presents the angles of these dendritic seg- 
ment arc-lengths for the alpha (a) and beta(b) cells. The 
last features considered in this work refer to the branch 
lengths and angles at the branch points, which are shown 
in Figures [7| and ODfor the alpha (a) and beta (b) cells. 
By "branch length" it is understood the total arc-length 
while moving from the branching point to the cell soma. 

Note that both branching angle densities (i.e. Fig- 
ures |S] and [TJ are similar for both alpha and beta cells. 
The length-related densities (i.e. Figures [S] and |HJ were 
obtained for alpha cells and then normalized with respect 
to their respective diameters (i.e. the largest distance 
between any two points of each cell) in such a way that 
they have the same average diameter as beta cells. Such 
a normalization was adopted so that the percolation only 
reflects the shape (and not size) of the cells. 

II. NEURONAL SYNTHESIS 

In order to generate the neuronal shapes, the proba- 
bilistic model of the neuronal geometry described in the 
previous section was statistically sampled by the Monte 
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FIG. 5: Distribution function of dendritic segments arc- 
lengths to the alpha (a) and beta (b) cells. 



Carlo approach as explained in the following. Initially, 
the soma of each cell was uniformly (Poisson distribution) 
distributed along an N x N matrix (associated to a dig- 
ital image) . The number of branches emerging from the 
soma was randomly chosen according to the respective 
density, being uniformly distributed along the somata, 
which are circular. For each cell, for each branch, the 
orientation of the emerging segment was drawn from the 
respective distribution. Straight segments are then in- 
corporated, piece-by-piece, into the growing process. The 
length and orientation of each of these segment pieces was 
sampled through Monte Carlo from the respective statis- 
tical model, therefore taking into account the previous 
angle and length. In order to allow all neuronal cells to 
grow in a 'simultaneous' fashion, a single segment piece 
is incorporated into each growing branch, for each neu- 
ronal cell, at a time ('round- robin' scheme). Every time a 
new branch was visited, the probability for new branch or 
growth termination was sampled, and the respective ac- 
tion taken. In case we have a branchp^|. the orientations 
of the two branching new segments were sampled from 
the respective distributions, and those branches were sub- 
sequently included in the 'round-robin' growth scheme. 
The growth of branches continued until one of the fol- 
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FIG. 6: Distribution function of dendritic segments arc-angles 
to the alpha (a) and beta (b) cells. 



lowing conditions is reached: (a) it is selected for inter- 
ruption; or (b) it reached 10 hierarchical stages. Fig- 
ure |5] illustrates morphologically-realistic neuronal net- 
works obtained by the growing process described above 
considering alpha (a) and beta (b) cells. 

III. PERCOLATION DYNAMICS 

During the simulated neuronal growth, a synaptic con- 
nection is implemented every time a growing dendrite 
overlaps any portion of the other current cells. So, as the 
cells develop in size and shape, they tend to form more 
connections. A group of connected cells is henceforth un- 
derstood as a cluster. A natural representation of such 
growing structures can be immediately obtained by using 
graphs whose nodes correspond to the neuronal cell soma 
and the edges correspond to the synaptic connections. 
While several topological and morphometrical properties 
of the evolving neuronal networks can be quantified, in 
this work attention is concentrated on the size S(t) of 
the cluster containing the maximum number of nodes — 
i.e. the dominating cluster — found at each time instant 
t (i.e. the growing stage). The sizes S(t) are calculated 
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FIG. 7: Distribution function of the angles at the branch point 
to the alpha (a) and beta (b) cells. 



from the graphs which are constructed as the networks 
evolve. The critical phenomenon of percolation is identi- 
fied by looking for an abrupt transition along S(t), which 
is related to the formation of the giant cluster Q . After 
this point, the growing neuronal structure is character- 
ized by the presence of such a giant community, which 
dominates the subsequent connectivity dynamics. 

A single model of cell was considered in each simula- 
tion in order to keep statistical variability low and allow 
a more precise identification of the percolation critical 
point (not a density as in traditional percolation the- 
ory, but a time instant during the neuronal outgrowth). 
The chosen neuronal model is "stamped" N times on 
the considered space (a rectangular window of 1000 by 
1000 elements) according to the uniform probability. A 
total of 500 realizations was performed for each consid- 
ered configuration, from which the average and standard 
deviation shown in the graphs were obtained. In order 
to avoid intense superposition between cells, cells were 
placed at least 5 pixels apart one another. 

Figure lTUT a-d) presents the evolution of the maximum 
cluster size considering growing densities of alpha cells, 
while Figure HUT e-hl presents analogous graphs consider- 
ing beta cells. As expected, the critical transition tends 
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FIG. 8: Distribution function of the branch lengths to the 
alpha (a) and beta (b) cells. 



to increase with the density of neurons, with markedly 
sharper transition being verified for the beta neuronal 
cells. More interestingly, it is also clear from the ob- 
tained results that the alpha cells implied longer times 
(e.g. higher cell densities) before percolation. Generally, 
percolation was often observed after 400 growth steps 
for beta cells, but only after 600 steps for alpha cells. 
This is exactly the opposite result than it would be ob- 
tained in case the cells were not size-normalized. Indeed, 
the fact that size-normalized beta cells tended to perco- 
late sooner than alpha cells provides a clear indication 
that such cells tend to have more intrincated morphol- 
ogy. As indicated in Figure which describes the cu- 
mulative two-variated distribution of dendritic segment 
angles, beta cells are characterized by higher dispersion 
of angles, implying the overall dendrites to become more 
disordered and spatially complex, which is in full agree- 
ment with the obtained percolation dynamics, i.e. more 
complex neuronal cells tend to percolate sooner than less 
complex cells with similar sizes. 
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FIG. 9: Examples of neuronal networks obtained by the adopted growth methodology considering the alpha (a) and beta (b) 
models. Different gray levels were used for neuronal cell representation in order to facilitate the visualization of the individual 
cells. 



IV. CONCLUDING REMARKS 

This article has reported on several new perspectives 
related to neuromorphic models and percolation induced 
by dendritic growth. First, we have shown how mor- 
phologically realistic neuronal networks can be simulated 
by using Monte-Carlo sampling of statistical models de- 
rived from a series of geometrical measurements of real 
neuronal cells. Second, we have investigated a new per- 
spective to percolation studies in which, instead of in- 
corporating new connections of fixed size between the 
involved elements, the percolation dynamics is defined 
by the progressive growth of dendrites/ axons, follow- 
ing biologically-realistic rules derived from experimental 
data. The obtained results indicate that the percola- 
tion in such evolving systems is also characterized by 
abrupt transitions of the dominating cluster size along 
the progression of the growth and connections. We have 
shown that distinct critical points are usually identified 
for growing dynamics of systems underlain by distinct 
neuronal morphologies, with beta cells reaching perco- 
lation sooner than alpha cells, a result that is related to 
the fact that the dendritic processes of beta cells are more 
intrincated and spatially complex. 

Such results establish and interesting connection be- 
tween the statistical geometrical features of the consid- 
ered cells and their potential for forming clusters among 
the neuronal milieu. Such perspectives and results are 



particularly interesting because the functional properties 
of neuronal networks are closely related to their connec- 
tivity (e.g. [HI [H HHIl). The perspectives for further 
investigations are many. For instance, it would be par- 
ticularly interesting to check how the consideration of 
more than one distinct statistical model of neuronal ge- 
ometry will affect the measured critical point. Another 
interesting possibility is to investigate, in the spirit of 
|l8j . to what an extent the critical point statistics can 
be used as a resource for classification of the morpholog- 
ical types of involved neuronal cells. A third promising 
future development is to quantify, through simulations, 
how the geometrical properties of the neuronal cells (e.g. 

by controlling the sizes of the neuronal clusters, ul- 
timately define the functional properties of the obtained 
structures 0. Still, it would be interesting to correlate 
several measurements from complex network research, es- 
pecially those related to the hierarchical structure of the 
networks [2(j, with the critical percolation time. 



Acknowledgments 

Luciano da F. Costa is grateful to FAPESP (processes 
99/12765-2) and CNPq (process 308231/03-1) for finan- 
cial support. The authors thank Luis Diambra, Marconi 
S. Barbosa and Gonzalo Travieso for reviewing and com- 
menting on this work. 




FIG. 10: Mean and standard deviation of the size of the largest cluster in 
cells) of alpha (a-d) and beta (e-h) cells. 



terms of increasing densities (300, 400, 500 and 600 
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